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Abstract — We derive a modular fluid-flow network conges- 
tion control model based on a law of fundamental nature in 
networks: the conservation of information. Network elements 
such as queues, users, and transmission channels and network 
performance indicators like sending/acknowledgement rates and 
delays are mathematically modelled by applying this law locally. 
Our contributions are twofold. First, we introduce a modular 
metamodel that is sufficiently generic to represent any network 
topology. The proposed model is composed of building blocks 
that implement mechanisms ignored by the existing ones, which 
can be recovered from exact reduction or approximation of 
this new model. Second, we provide a novel classification of 
previously proposed models in the literature and show that they 
are often not capable of capturing the transient behavior of the 
network precisely. Numerical results obtained from packet-level 
simulations demonstrate the accuracy of the proposed model. 

Index Terms — congestion control modeling; fluid-flow model; 
queuing model; self-clocking 

I. Introductory discussions 

A. The congestion control problem 

The congestion problem [[TJ , J2) is inherent to communica- 
tion networks where capacity of supporting infrastructure that 
relays information is small compared to user demand. Conges- 
tion is responsible for delay and data loss, which compromise 
the efficiency of the overall network. Controlling congestion is 
hence an important problem for which several algorithms have 
been developed. They mainly rely on the concept of congestion 
window^, which is adapted according to a congestion measure. 
According to the type of congestion measure, two classes 
of congestion control algorithms may be identified 0. El- 
The first and oldest class is loss-based, meaning that the 
congestion measure is the packet-loss information. This class 
is easy to implement but leads to a quite rough control 
since the protocol detects the network congestion only after 
provoking it. In order to control congestion more smoothly and 
prevent data loss, delay-based algorithms, using for instance 
the Round-Trip Time (RTT) information as the congestion 
measure, may be considered instead. They are however more 
difficult to implement due to the possible unavailability of 
certain necessary measures, such as queuing delays. 

The main difficulty in congestion control lies in the fact 
that, basically, the hosts ignore almost everything about the 
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network: the routes and their capacity, the numbers of routers 
(hops), the number of users, etc. Hence, protocol designers 
face the problem of controlling a very large and complex 
system with actually very little information. 

When designing a protocol, stability of the network is 
certainly the most important constraint. Performance criteria 
can be additionally considered in order to optimize the network 
behavior. For instance, we may want to use all available band- 
width {efficiency), share it equally between users (fairness), 
and/or be tolerant with respect to unregulated traffic and other 
protocols (cross-traffic adaptation). 

In order to observe/predict the network behavior and va- 
lidate a protocol, simulations and experiments must usually 
be conducted. NS-2 is a widely accepted open-source event- 
based simulator dedicated to this purpose. However, as any 
other simulator, it does not permit to analyze the stability 
of a network theoretically. Hence constructing mathematical 
models for networks may play an important role in network 
analysis and protocol design since they potentially allow for 
a theoretical analysis and an equation-based design of new 
protocols. 

B. Models, approximations and accuracy 

Modeling is now ubiquitous. The key idea is to start from 
a system and arrive at an abstract representation of it, such as 
one given in terms of a set of mathematical equations. It is not 
always necessary that a model represents all characteristics of 
a system but only a subset of interest: e.g. a molecular-level 
model can be irrelevant to portray a river. This gave rise to 
fluid-mechanics which, although being an idealization of the 
reality, yields very accurate predictions. A similar idealization 
has been shown to be very useful for the congestion control 
problem through the consideration of fluid-flow models [4]. 

C. Metamodels and network models 

A paragon of network modeling is undoubtedly used in the 
field of electrical engineering. It is easy to identify the reasons 
for the success of the theoretical framework: 

1) Only two universal concepts: current and voltage, gov- 
erned by simple laws (Kirchhoff's laws) 

2) Local description of the elements in terms of these 
variables and additional local concepts (e.g. resistance, 
etc.) 
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3) Easy transcription of the electrical network into a topo- 
logically identical diagram, and vice-versa. 

4) New models corresponding to new devices may be freely 
added without compromising the existing ones. 

5) Model predictions fit very well to reality. 

6) Systematic way of analysis by hand calculations or 
simulators. 

A very important feature is that the principles of modeling an 
electrical network is independent of its topology and elements. 
This is achieved thanks to the structure of the paradigm that 
we refer from now on as a metamodel, which is a model that 
consists of a set of frames, rules, constraints, submodels and 
theories applicable and useful for modeling a predefined class 
of problems. In the case of electrical networks, the metamodel 
consists of the concepts of current and voltage, the Kirchhoff 's 
laws and the local models of electrical elements (resistor, ca- 
pacitor, transistor, etc.) as well as all the related mathematical 
tools. Since networks (like communication networks, electrical 
networks, transportation networks and even social networks) 
consist of interconnections of several elements, it turns out 
that metamodels are then very suitable for describing them 
since they also consist of interconnection of concepts, rules 
and submodels. Hence metamodels provide, in essence, an 
elegant scalable and modular way for modeling networks. 

D. Motivations and contributions 

The main motivation of this work is to give a clear picture 
of congestion control modeling problem through derivation 
of a metamodel having solid mathematical foundations. We 
introduce a modular metamodel that would lead to an inter- 
esting step forward towards a generic way of providing models 
for communication networks. This metamodel should then 
satisfy the additional constraints on independence of network 
topology (scalability) and elements (modularity). It should 
also provide accurate predictions along with simple graph 
representation. An underlying difficulty is the presence of sev- 
eral phenomena at different levels: decision to send a packet, 
transmission of packets on transmission channels, storage of 
packets in queues and time-varying waiting-time (queuing 
delays), congestion window size adaptation, etc. Finding a 
unified way for representing all these critical phenomena is 
challenging. 

The proposed metamodel is based on a single concept of 
information conservation, from which models for the network 
constituents (i.e. transmission channels, queues and users) 
are obtained. This allows to derive new models for network 
elements, obtain mathematical proofs for unproved/claimed 
existing ones, and invalidate some of them. All important 
variables of the network (sending rates, ACK rates, queue size, 
etc) are described by explicit formulas, hence computable. 
Using the proposed metamodel, describing a given topology 
is immediate, and performed by simply plugging the models 
together, so as the actual topology is transcribed into a graph 
having the different network elements located on the edges, as 
in electrical engineering. It is also proved that existing sending 
rate models are either approximation of the proposed sending 
rate model, or even exact provided that the network topology 



satisfies certain structural conditions. Simulations and com- 
parisons with existing works tend to suggest the relevance, 
reliability and accuracy of the proposed metamodel. A non- 
exhaustive summary of related works on congestion control 
modeling is finally made in order to compare congestion 
control models according to important properties and criteria^. 

The outline of the paper is as follows: Section [TT] intro- 
duces the particular network graph representation considered 
in the paper. Using continuous-time models, such as fluid- 
flow models, to describe networks is justified and concepts 
of universal clock, local discrete-time network element clock, 
and clock-coupling are defined in section [HI] Section [IV] 
presents the conservation law of information and we derive 
the transmission channel model in section [V] In section [VT] 
we develop the model for the FIFO buffer network element. 
The user model is given in section IVIII and we summarize the 
obtained results in a compact form in section IVIIII In section 
IIX1 we consider a network with single buffer/multiple-user 
topology to implement the proposed model. We validate our 
model in section [X] and related work is given in section IXll 
Section IXTTl concludes the paper. 

II. Networks and Graphs 

It is convenient to introduce here the particular network 
graph representation considered in the paper. It is different 
from the traditional ones 0, @, Q, |H] since it places all 
network elements on graph edges, leaving nodes with the role 
of connecting points, as in electrical circuits. Four types of 
nodes are distinguished: the input nodes b J and output 
nodes uf, for user m and buffer bj respectively. The 
superscripts have to be understood as a temporal order of 
reaction or causality: the data come (-) then leave (+). We 
will denote any edge E of the graph by (x, y) where x and y 
are the input and output nodes respectively. Moreover, given 
any edge E, the input and output nodes are given by j3(E) 
and e(E) respectively. 

According to these definitions, a queue edge is always 
denoted by (b~,bf), a user edge by (u^,it+) and a trans- 
mission edge by (bf,uj), (u~l,bj) or {bf,b^), i ^ k. This 
is illustrated in Fig. Q] We call a circuit C; the communication 
path of user Ui, that is the path connecting its output u~l to its 
input u~, i.e. Ci — (u^,u~). Note that in complex networks 
there exist several possible paths but only one of them, the 
one used for communication, is a circuit. In Fig. Q] the only 
possible circuit is given by C = {u + , b~ , 6 + , u~). 

III. Fluid-flow idealization 

In an asynchronous network like the Internet, each el- 
ement can be considered to have its own local discrete- 
time clock Tj C K.+ where Tj is countable, governing 
the rhythm of protocol decisions and packets transmission. 
In congestion control, the clocks beat with the rhythms of 
acknowledgment reception rates, which are influenced in turn 
by network congestion; this is referred to as ACK-clockin^ 

- The authors are aware of the fact that all the models are certainly not 
listed. This is however a first attempt and any suggestion from the reviewers 
to improve/correct this comparison is welcome. 3 The term self-clocking is 
also used but is less explicit. 
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(6", 6+) 

b~ r — b+ 

(u + ,b~) ! h \ \ i2 \ (b+,u-) 

Fig. 1. Example of graph with 4 edges: one user edge u = (u~ ,u + ), 
one buffer edge b = (b~ ,b + ) and two transmission edge l\ = (u+,b~), 
£ 2 = {b+,u-} 

When several sources send data through the same buffer/path, 
a flow-coupling takes place leading then to clock-coupling. 
This clock-coupling arises at a very large-scale and distant 
sources having their clocks coupled cannot be considered to 
have independent behaviors. As a consequence, the sending 
rates and acknowledgment rates are hence intimately inter- 
dependent. Modeling this clock-coupling and the underlying 
phenomena is of incredible complexity since the number 
of clocks and their interactions grow very quickly with the 
network complexity, leading then to a very complicated struc- 
ture for the interrelated local clocks Tj, see for example |9] 
Equations (3.7)]. 

An idea to resolve this complex time-structural problem 
relies on the definition of a universal clock T u dictating a 
common time to the entire network. This leads us to the 
following fact: 

Fact 1. There exists an ideal universal clock T u embedding 
any local clock T 4 , i.e. (\J. Tj) C T". 

A natural universal clock is given by |J . Tj and is a discrete- 
time clock. It however does not simplify too much the mod- 
eling problem since it is difficult to write recurrence relations 
for general network topologies [9 |. Deriving a metamodel 
achieving scalability is then unlikely using such a universal 
clock. A however less natural universal clock T" assimilated 
to a clock running over positive real numbers continuously, i.e. 
T" = K + , is much more promising. This particular universal 
clock indeed dramatically simplifies the modeling problem, 
and this motivates its consideration in this paper. Using such 
a time-scale, a metamodel can be obtained, resulting then 
in a scalable solution in which the network asynchrony is 
captured through appropriate expansions and compressions 
of the time-space. Furthermore it enables the use of well- 
established mathematical tools: real function analysis, integra- 
tion theory, dynamical systems, delay-differential equations, 
etc. A conclusion is that continuous-time models may be 
used to describe networks iflOl . ||4], ifTTI . Ifl2l . j6): these are 
generally referred to as fluid-flow models, emphasizing the 
connection with continuum mechanics and more specifically 
with fluid-mechanics. 

Within this framework, it is possible to provide a proper 
definition for data flows. 

Definition 2. Let us consider a point x in an edge E of 
the network and denote the number of packets having passed 



through point x between £o S and £ € T", £ > £rj, by 
N x (t,to), Then, the flow of data passing through point x is 
defined as a function cj> : E x TP 1 — > R + verifying: 

N x (t,t Q ) = f (f>{x,s)ds (1) 

J to 

where the integral is a standard one, e.g. the Lebesgue integral. 

Flows are defined in such a way rather than being the 
derivative of the number of packets, since the number of 
packets is non-differentiable, i.e. flows may contain dirac 
pulses, steps and so on. It is also interesting to note that 
since the the universal clock embeds all the local clocks, it 
is possible to recover discrete-time asynchronous models such 
as the one in J9) by setting flows to be trains of dirac pulses 
on U, 

Using the notation defined in Section [TU we can build the 
flow vectors <fi(x,t) using the 'col' operator: 

<t>(x,t) = col [</> k (x,t)] (2) 

k= 1 

where x in any input and output node of the network elements, 
i.e. x can be any u~ , uf, b~ , bf. The quantity r](x) denotes 
the number of flows passing through node x. The concept 
of flows of data is hence very close to those of current in 
electrical engineering and flow of a liquid in fluid-mechanics. 

IV. Conservation law of information 

The core of the metamodel is the conservation law of 
information stated in this section. This law allows to improve 
the characterization of the elements of the network by notably 
clarifying their input/output relationship, enabling then a 
modular formalism. This conservation law follows from the 
remark that the quantity of information^ is preserved in a 
communication network: the data can either be in transit, 
lost or received. Assuming lossless networks, it is possible to 
determine the total number of packets in transit in any edge, 
simply by counting the number of entering packets according 
to a simple rule. When applied to a specific element, this 
law allows to characterize the fact that the information is 
preserved from the input to the output. 

Law of Conservation of Information: 

Given any edge E of a network, then for all ( 6 T", there 
exists a time £rj(£) € T", to(t) < t such that 

Ps(t) ■■= L4>{o,t)de 

= Jl (t) cf>(0(E),s)ds (3) 

= Np (E) (t,t {t)) 

The integration over E is an abstract integral which has 
to be understood as a flow spatial integration from (3(E) 
to e(E), that is, the number of packets Pe(*) m the edge 
E = {(3{E),e{E)) at time t. 

The above result stated in quite abstract terms just says that 
the number of packets in transit in an edge at a certain time t 
can be determined by counting the number of entering packets 
(i.e. integrating the input flow) over the interval [£n (£),£], the 

4 expressed in bits or packets. 
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Constant propagation delay T 




Fig. 2. Transmission channel block 

lower bound to of the interval depending on the considered 
element, i.e. transmission channel, queue or user. A simple 
application of the law is given in the next section discussing 
transmission channels models. 

The first benefit of this law is to show that we can inter- 
changeably use a spatial or a temporal integral to calculate 
the quantity of information (number of packets) in an edge. 
The temporal integral formulation is very convenient to work 
with since it requires the knowledge of the input flow only, 
rather than the flow value on the entire edge for the spatial 
integral formulation. This hence allows to discretize the space 
dimension by only considering flows at the nodes, simplifying 
then the network representation and the modeling problem. 

The second benefit lies in the fact that the temporal integral 
can be utilized to yield explicit solutions for the output flows 
of the different network elements. A very general result is 
given below: 

Proposition 3. The input flow <p(/3(E), •) and the output flow 
(j>(e(E), •) of edge E verify 

4>(e{E),t) =t (t)'4>(P(E),t (t)) (4) 

where we assume that to(t) is absolutely continuous and 
to(t)' is the upper-right Dini derivative of to at t, i.e. 
t' (t) = limsup^o hr 1 (to(t + h) - t (t)). 

Proof: Since Np/E)(t,to(t)) is the current number of 
packets on edge E at time t, then differentiation with respect 
to time provides the balance equation 

[N 0(E )(t, t (t))}' = 4>{P(E),t) - to(t)'0(P(E),t o (t)). 

Note also that a second valid balance equation is given by 

{N p{E) (t,to(t))Y = WE),t) - <t>(e{E),t). 

Identifying the right-hand side yields the result. ■ 
The proposition given above plays a crucial role in the 
metamodeling problem since it provides an explicit formula 
of the output flows. This output flow verifies the conservation 
of information from the input to the output of the edge E. 
By integrating the input and output over [0, oo) the very 
same value is obtained. This emphasizes that the output flow 
is defined in such a way that, as desired, it respects the 
natural property of conservation of information. Proposition [3] 
is used repeatedly in the paper in order to provide accurate and 
explicit models for transmission channels, queues and users. 
Applications are given in Sections [V] IVI-CI and IVII-DI 

V. Transmission Channel Model With Constant 
Propagation Delay 

The first element-model is derived in this section, namely 
the model for lossless transmission channels with constant 



delay. They exactly behave as transmission lines and a delay- 
based formulation is provided. The derivation is rather straight- 
forward but it is a good example of application of Proposi- 
tion [3] 

Result 4. Given a lossless transmission channel correspond- 
ing to edge E and having constant propagation delay T > 0, 
the output flow is given by 

<t>(e(E),t) = <f>(l3(E),t-T). (5) 

The corresponding module is depicted in Fig. [2] 

Proof: According to the conservation law (|3), the number 
of packets in transit Pe (t) in the edge E at time teT" obeys 



<p(/3(E),s)ds 



(6) 



to(t) 



where to (t) = t — T since the propagation delay is constant. 
A packet sent at time t — T will indeed be, at time t, still 
in the edge and about to leave. The result follows then from 
Proposition [3] ■ 
We are now in a position to introduce transmission channel 
operators defining part of the network topology. 

Definition 5. The flow vectors 4>{u~ ,t) = col^ [4>{uJ , t)], 
<j>(b~,t) — co\i[(j)(b~ ,t)}, (j)(u + ,t) — coli[(f>(uf ,t)] and 
<f>(b + ,i) = coU[(j)(bf , t)} are related by transmission channel 
operators lZ m , • G {ub, bu, bb} as 



where 



K 



(u-,t) 
(b-,t) 







K 



Hub 
T^bb 



cj>(u+,t) 
0(6+, i) 



and V 



VS(t) 



(7) 







The matrices lZ m , • G {ub, bu, bb} correspond to routing 
matrices between output and input nodes. For instance, lZ u b 
maps flows at user output nodes to flows at buffer input nodes. 
These matrices essentially consist of constant delay operators 
with delays corresponding to transmission channels. The full- 
rank input matrix T>b drives the vector of cross-traffic flows 
S(t) to buffer input nodes. 

VI. FIFO Buffer Model 

This section is devoted to the very important buffer element 
which temporarily stores incoming information before pro- 
cessing it. First, the standard fluid model for queues is recalled 
ifTTI . (3), and a complete delay-map is characterized. In 
order assign each input flow to its corresponding output flow 
and solve the output flow separation problem iTPJl . iTPfl . lfT31 . 
ESI, IfTTI . ifTHl . the conservation law is then applied on the 
standard queue model, in a similar way as for transmission 
channels. The output separation problem allows us to focus 
on the accurate description of queues, which captures both 
the FIFO behavior and the form of output flows. Some extra 
discussions and interpretations of the results are also provided. 
Finally, a comparison to an existing model for output flows is 
carried out and concludes in favor of the proposed one. 
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A. Queue model 

Routers have queues to store incoming packets temporarily. 
The following integrator model ifTD can be proved to be a 
limit model of a M/M/l queue when the packet size and thus 
the processing time tend to J5], |[T9l . 

Definition 6. The queue dynamics of buffer bi is governed by 
the model 

*(*) = (8) 

3=1 

where the aggregate output rate is defined as 



Ci if Ci(t) 

J2j 4>j (K ' *) otherwise. 



(9) 



Above, gj, Ci and (f>j(b~,t) represent the queue size, the 
maximal output capacity and the flow of type j at the input 
respectively. The condition Ci(t) is given by 



Ci(t) 



h(t) > o] v 



fa ( b * > *) > Q 



(10) 



The corresponding queuing delay can be easily deduced using 
the relation Tift) — qi(t)/ci. 

The above model can also be refined to capture additional 
features such as finite maximal queue length, flow priorities, 
and multiple output capacities. These extensions are omitted 
here since they are straightforward. It is important to stress 
that this model is incomplete and useless in this form. First, 
the output flow is given in aggregate form. This prevents 
the modeling of the appropriate routing of each output flow. 
Second, it does not capture the queue FIFO behavior. Finally, 
the model does not assign a specific queueing time to each 
input flow. In the following subsection, the conservation law 
(0 is used in order to confer the FIFO property to the model 
and separate the aggregate output flow into distinct flows. 

B. Forward and Backward Queuing Delays 

The maps defined in this section are very useful for obtain- 
ing a closed formula for the buffer output flows in Section 
|VFC] and for RTT in Section IVlFCl 

Let us consider first the buffer model dg) with queueing 
delay Tj(i). Assume that the time instants at which packets 
enter the queue are chosen as reference times. We may then 
be interested in predicting the packet output time. This leads 
to the following definition: 

Definition 7 (Forward delay operator). The forward delay 
operator fi : T u — » T" corresponding to buffer bi mapping, 
at a flow level, any input-time t to the output-time fi(t) is 
defined as 

fi(t):=t + n(t) (11) 

where tj (t) is the queuing delay of buffer 

It is easy to see that output-time can be readily computed 
from the knowledge of input- time and queuing delay value. 
If however, we would like to set the reference time to be the 



output time, we may ask the question whether it is possible 
or not to retrieve the input-time from it. This is equivalent to 
asking the question of invertibility of the map f. 

Result 8 ([ 17 1). The map fi is invertible if and only if the input 
flow of the corresponding buffer is positive almost everywhere. 
▲ 

Hence, provided that there is a nonzero input flow to the 
buffer, the input-time corresponding to a given output-time is 
well-defined and can be obtained using the backward delay 
operator. 

Definition 9. The backward delay operator gi : T" — > T" 
corresponding to buffer bi mapping, at a flow level, any output- 
time t to the input-time gi(t) is defined as gi := ff under 
the assumption of Result^ 

We also have the following useful results: 

Result 10 ( 1171 ). The functions fi and gi := ff obey: 



m 



m = 



[n(9i(t))}' 



<t>k(b l ,gi(t)) 



Ci 



if Ci(gt(t)) 
otherwise 

if Ci( gi (t)) 
otherwise 

if Ci (<?»(*)) 
otherwise 



where f'(f) stands for the the upper right Dini derivative of 
f(t), i.e. f'(t) = limsup, 40 h- 1 (f(t + h)~ /(*)) . A 

The following technical result allows to simplify the condi- 
tions involved in hybrid models: 

Result 11. The equivalence C(gi(t)) C(f) holds. 

Proof: Proof of =>: If the buffer is congested at time 
gi(t) then the buffer will also be congested at time t since the 
data entered at time gi(t) leave at time t. 
Proof of <^: Conversely, if there is any data to leave at time 
t, they must have entered in the queue in the past, i.e. at time 
gi(t). Equivalence is proved. ■ 



C. FIFO Buffer Output Flow Separation 

In this section, we use the results given above and the 
conservation law in order to improve the buffer modeling by 
adding the FIFO characteristics and splitting the aggregate 
output flows into a sum of distinct ones. Without further 
consideration on the queue type, there exists an infinite number 
of ways to separate the aggregate output flow directly from 
the queuing model of Definition [6] When a FIFO queue 
(i.e. order preserving) is considered, it turns out that the 
output flow separation problem is easily solvable. The FIFO 
characterization and output flow separation problems have 
been fully solved in fP71 . In this section, we will simply recall 
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and explain these results, and connect them to the conservation 
law ®. 

Result 12 ( 11171 ). Let us consider the queueing model (0 
which we assume to represent a FIFO queue. The output flow 
corresponding to the input flow 4>i(b^,t), I < rj{b^) is given 
by 



Forward queueing delay r 



Cjfajbj ,9i(t)) 

Ef i] <t>j(K,9i(t)) 



if Ci(t) 
otherwise 



(13) 



Proof: A proof is available in 0171 . A more direct one 
based on Proposition [3] is given here. Noting that for buffer 
hi, we have to(t) = Qi(t), then using Proposition (0 and the 
formulas of Result [10] we get equation ( fT3l with the difference 
that the condition is C(gi(t)). However, from Result [TT] the 
condition C(gi(t)) is equivalent to C(t) and the result follows. 

■ 

The same model has been also proposed in fl3l . Ifl6l but 
claimed without any proof. We have shown above that this 
model is an immediate consequence of the information con- 
servation law and gives, for the first time, a theoretical proof 
for it. This considerably strengthens the trust we may have in 
this model. A comparison with packet-level simulations in the 
next subsection tends to show its exactness. 

This model also deserves interpretation: formula (Q~3) says 
that output flows consist of scaling and shifting of the input 
flows. The delay accounts for high flow viscosity and captures 
the queue FIFO behavior, at a flow level, while the nonlinear 
ratio expresses the flow coupling at the core of the flow and 
clock-coupling phenomena, see section [Till since each output 
flow depends on the corresponding input flow and all the 
other ones as well. A change in a single flow will affect 
all the output flows. This model also tells that the output 
flow corresponding to the input flow <fit(b~ ,t) is expressed 
as a (delayed) ratio of the input flow (pi(b~,t) to the total 
input flow that entered the buffer at the same time. Hence, the 
output flows are proportional to relative flows modeling the 
'chance' of having a packet of certain type served at time t. 
This 'chance' is then scaled-up by the maximal output capacity 
to utilize the available capacity. 

Operators representing buffers can now be introduced: 

Definition 13. The buffer operator Bi with r\(bi) input flows 
is defined as 



B, 



<Kb+,t) 



(14) 



where the output flows and the buffer state are governed by 
dSJ) and ( 1751 ). Using these operators, we can build a matrix of 
operators B connecting the (p(b~ ,t)'s to the (f>(b + ,t)'s as 



cl>{b + ,t) = B4>{b-,t) 



(15) 



where B = diag^ {Bi}. 



Queue (b~,b+) 



- 0(b+,t) 
defined in J 1 3t 



Backward queueing delay r(g) 
Fig. 3. Queue/Buffer block 

D. Comparison with Another Model 

Two main models for buffer output flows have been reported 
in the literature on fluid-flow models: the flow-based model 
iTOl . lTl6l . 11171 described in this paper and the pseudo-queue- 
based one Oil, 03), Q8) given by 



E*«?(*) 





if ft(t)^0 
otherwise 



(16) 



where qf(t) is the number of packets of type £ in queue i and 
Ci is the maximal output capacity of queue i. 

Until now, these models have not been confronted to each 
others. In the following, they will be theoretically and exper- 
imentally compared, and it will be shown that the flow-based 
model is the only model that faithfully characterizes the actual 
output flows, validating then the proposed conservation-law- 
based paradigm. 

1) Theoretical argumentation: First, the above model as- 
sumes that the 'chance' of having a packet of type I at the 
output at time t is q\{t)j g| (t). Model ([ToT) then makes no 
difference in picking a packet in the middle, at the end or at 
the beginning of the queue since only the number of packets 
matters. It is thus unable to capture the FIFO characteristic 
of the queue since swapping packets in the queue does not 
modify the output flow. In contrast, the proposed model does 
capture this characteristic through the delay dynamical model 
and the delayed nonlinear input-output relationship involving 
flows directly: relative variations of the input flows are passed 
to the output flows after some queueing delay. Note however 
that both models coincide at equilibrium. 

Second, since output flows in model dT6l are computed 
from the integration of input flows, it turns out that the 
map from input flows to output flows is a nonlinear low- 
pass filter with 'bandwidth' equal to l/q(t) when q(t) > 0. 
High frequencies in the input flows are hence filtered out, 
making the existing model inaccurate for high frequency flows 
(fast transient), especially when the queue size is large. Note 
that the actual buffer behavior does not have any filtering 
effect, it just behaves as an ordered tank. On the other hand, 
the proposed model does not filter out any frequency band 
due to its feedthrough structure. It however has a distortion 
effect on the input flows due to the nonlinear structure and 
the dynamically changing delay, i.e. change of frequency and 
amplitude. As a result, the queue-based model does not satisfy 
any conservation law since low-pass filters dissipate energy all 
over the frequency band, resulting in information loss at the 
model level. This is in total contradiction with the actual queue 
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behavior that just stores information and does not dissipate 
anything. Note however that the flow-based model intrinsically 
satisfies the conservation since the model is derived from it. 

Last, the proposed model incorporates naturally the queuing 
delay in the expression, while for model (fTol it is unclear 
which delay to consider, since the order of arrival of data is 
not tracked. 

2) Case Analysis: To compare the models, let us consider 
two input flows given by 



(i 
(i 



/3)c(l 
/3)c(l 



Sq(wi)) /2, 
Sq(wi)) /2 



(17) 



where c > and u> > are the buffer output capacity 
and the oscillation of flows. The term (3 > is a tunable 
parameter related to the amplitude of the inputs flows and 
the function Sq(utf) := sign(sin(W)) is a square function of 
period T := 2tt/ui. Since the flows are in phase opposition, 
they lead to an alternation of packet types in the queue while 
packet populations remain roughly close to each other at any 
time. Therefore, the output flows should reflect the actual 
content of the queues and the model should be able to keep 
track of the order of arrival of packets in the queue. 
The queue-based model (IT6b predicts the queues 



1//3 



9< (t) = 0(t)fc(O) + / 0(t-a)&(«)ds, * = 1,2 (18) 
?(0) 



with 0(t) = 



from which it is quite difficult 



to foresee the shape of the output flows. We can however note 
that the filtering effect of the convolution operator with kernel 
9 will deform the input flows, making then the output flows 
not square anymore. 

When the proposed flow-based model is considered, it is 
enough to compute the forward and backward delays, and 
apply the formula for output flows: 



r(t) 
/(*) 

m 

Mb + ,t) 



(1 



1 



r(0), 

;(*-t(0)), 



Mb + ,t) = - 



1 + 

c(l + Sq(w5(t)))/2, 
c 

2 

c 

2 



(19) 




In this case, the predicted output flows have the same shape as 
the input flows but with a different frequency and amplitude. 
The proposed output flow model then exactly captures the 
content of the queue, that is the alternation of blocks of size 
N(T) of type 1 and type 2. It is also immediate to see that 
the number of packets N(T) = JVf(T), i = 1,2 received by 
the buffer over one period T is given by 

cT 



N(T) := / Mb + ,s)ds = 



7t(/3 + 1)c 



(20) 



By virtue of the information conservation law, the same 
number of packets is retrieved at the output over the period 
(j3 + 1)T, enlarged due to the limiting output capacity c. It 



is quite convincing that the flow-based model yields a much 
more coherent picture for this example. 

3) Simulation: For simulation, we choose lo = 2n (i.e. 
T = Is), /3 = 1, c = lOOMb/s and r(0) = 0. The output 
flows obtained from the different models are depicted in 
Fig. HJ where we observe the behaviors predicted by the above 
calculations. Notably, the output flows predicted by model 
( [Tol l tend to slow-down (low pass-filtering effect), decrease 
along time and seem to both tend to c/2, which is basically 
unrepresentative of the actual content and output of the queue. 
We can also notice the decrease of the bandwidth for the 
queue-based model as long as the queue grows in size. 

For comparison with NS-2, which deals with packets rather 
than flows, we integrate (up to an additional constant) the 
output flows predicted by each model to obtain a number 
of packets so that the comparison with NS-2 makes sense. 
The results are depicted in Fig. where we can see that 
the proposed model yields exactly the same results as NS- 
2 while the queue-based model is unable to track the stair-like 
curve returned by NS-2. It is also important to stress that the 
considered scenario is quite convenient for the queue-based 
model since the input flows contain mostly constant parts (low 
frequency parts). A very fluctuating input flow would be very 
penalizing and would make the low-pass filtering effect of the 
queue-based model even more apparent. 





Fig. 4. Comparison of the output Fig. 5. Comparison with NS-2 simu- 
fiows predicted by model 1131 (top) lations of the total number of packets 
and model )16t (bottom). Plain: output counted at the inputs of the buffers 
flow </>i(ft + ,t), dashed: input flow described by the proposed flow based 
01 (b~ , t) model t8l- )131 and the queue based 

model (U-GD- 



VII. Complete user model - Window control 

The derivation of the user model is, partially, still an open 
question and a complete solution, based on the conservation 
law (f3]), is proposed in this section. We assume that the 
implemented congestion control protocol admits a fluid-flow 
approximation, for instance interpolating the discrete-time 
trajectories of the real protocol [11], [20]. The conservation 
law is then applied over the circuit used by a user to obtain 
the so-called ACK-clocking model l2fl . E3 relating flight- 
size, user sending rate and RTT together. Several expressions 
for RTT are provided according to the considered reference 
time, similarly as for buffers. These results are finally merged 
together in order to clarify the connection between the flight- 
size and the user sending rate. The last step concerns the 
derivation of formulas relating the above variables to user 
congestion window size. 



■s 



A. Protocol model 

Here we assume that the congestion control algorithm 
can be represented as a continuous-time (hybrid) dynamical 
system. That is we have the following fact: 

Fact 14. There exist bounded functionals Vi, Wj and Ui such 
that the trajectories (zi(t),Wi(t)) of the following continuous- 
time model defined over T u 

Wi [t) = Wi(zi(t),tii{t)) (21) 
(j>i{uj ,t) = Ui(wi(t),4>i(u^,t)) 

match the trajectories of the asynchronous protocol ( defined 
on TjJ at points in T" Pi TV Above, Zi, Ui, 4>i(u^,-) and 
4>i(uf , •) are the state of the protocol, the congestion measure, 
the acknowledgment flow rate and the user sending flow 
respectively. The window size, Wi, is considered as the number 
of outstanding packets to track and supposed to be (weakly) 
differentiable. 

A procedure to solve the above interpolation problem has 
been first proposed in iTTTI for TCP and has been adapted to 
FAST-TCP in flU Appendix CJ. 

B. The ACK-Clocking model 

The ACK-clocking model ||2TI . l22ll is a very important con- 
sequence of the conservation-law. It characterizes the flight- 
siz<^ Fi{Ci,t) := Pd(t) °f the user Ui at any time t G T" 
over a circuit Ci = (uf,u^). The importance of the ACK- 
clocking model lies in the semantic; it adds to the model 
by relating RTT, flow and flight-size togethei[§ With this 
result, it is possible to incorporate a number of important user 
properties and mechanisms into the corresponding model. 

Result 15 (FS-ACK-Clocking). The (FS)ACK clocking model 
is given by 

Fi{Ci, t + RTTi{t}) = [ <P(9,t + RTT t {t})d6 

t+RTTi{*} 

= J <t>i(uf,s)ds 

(22) 

where RTTi{F\ is the RTT of a packet sent at time t by user 
Ui over the circuit Ci. k 

Proof: Since flight-size is a number of packets in a circuit, 
it can be cast as a spatial integration over the corresponding 
circuit. Thus, according to the conservation-law, it is possible 
to convert the spatial integration into a temporal one provided 
that we can determine the integration bounds. To obtain them, 
we use the notion of RTT and suppose that a data sent by 
user Ui in the circuit Ci at time t has a round-trip-time given 
by RTTi{£}. This means that the packets sent between t and 
t + RTT i{t} are unacknowledged at t + RTT t {t} and thus 

5 The number of outstanding packets. 6 Note however that in 1221 the 
window size is considered instead of the flight-size, which is rather different. 
To make the distinction, the window-based ACK clocking model is denoted 
by W-ACK while the fligh-size-based one by FS-ACK. Equivalence holds 
when some conditions, such as 'fiight-size>window size', are met. We will 
come back on this in Section IVlI-EI 



still in the circuit. Hence, the corresponding temporal integral 
has bounds t and t + RTT.^i}. ■ 

C. Round-Trip-Time Models 

The RTT consists of sum of a constant and a time-varying 
part, namely the propagation delays and the queuing delays, 
which have been characterized in Sections [V] and [VI] respec- 
tively. By combining these results together, it is immediate 
to obtain a RTT based forward model on a forward circuit 
operator, which considers the packet input time as a reference. 

To properly define it, let us consider a circuit C with N 
queues, indexed from 1 to N. The indices and N+l are used 
to denote the input and output nodes of the circuit respectively. 
Given a packet input time t, the corresponding packet output 
time tc and RTT obey the following formulas based on the 
forward circuit operator ,^c- 

Definition 16. The forward circuit operator &c '■ T" — > T" 
of circuit C assigned to any packet at input time t, an output- 
time tc(t) is given by 

tc{t) = ^c(t), 

&c = V^o/jjoV^o/^o... (23) 

°V^g o f 2 o /l Vq j 

where Vij is the constant delay operator with delay corre- 
sponding to the propagation delay between i and j and o the 
composition operator. The corresponding RTT expression is 
then given by 

RTT{t} =t c (t)-t= {&c ~ Id) (t) (24) 

where Id is the identity operator. 

Formula (l23l actually represents the alternation between 
constant delay operators corresponding to transmission chan- 
nels delay (the Vj's) and the queuing delays corresponding 
to queues (the fi i). Example [18] illustrates this formula on 
the topology depicted in Fig. Q] The same formulas, albeit 
expressed in different ways, have been also obtained in (9] 
Section 3.3.5] and 1221 equations (7d-7f)]. 

Although immediate to obtain, these expressions suffer from 
several drawbacks. First, operator &c lS clearly noncausal 
since it requires the knowledge of future information. Second, 
as pointed out in Section [VI-Bl the most convenient reference- 
time to use is the output time. In the user modeling problem, it 
coincides with the reception time of acknowledgments and the 
moment when the user receives the RTT information. Hence, 
it seems to be more convenient to consider a backward circuit 
operator based on the backward delay operator of Section 
IVI-BI The existence of such an operator is immediately 
inferred from the existence of the backward delay operator. 

Definition 17. The backward circuit operator 3§c ■ T U T U of 
circuit C assigned to any packet at output time tc, an input- 
time t(t c ) is given by 

t(t c ) = &c(tc), (25) 
S$ c = V ,i o gi o Vi, 2 ° 92° V 2 , 3 o ■ ■ ■ 
°9N-1 ^N-l,N 9N V w ,w+i- 
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Moreover, the corresponding RTT expression is given by 

RTT{t} =tc- t{t c ) = (Id - &c)(tc). (26) 

It is important to stress that the RTT formula given above 
looks noncausal since the RTT of a packet sent at time t 
is defined in terms of time tc > t. This is however not a 
problem since the RTT information is only available, and then 
used, by the user at time tc, when the ACK packet is actually 
received. What is important is the causality of the operator S$c 
in order to ensure computability of the RTT at any time. This 
emphasizes once again the relevance of considering output 
times as references. The example below illustrates the above 
discussions: 

Example 18. In the single-user/single-buffer case, the above 
expressions reduce to 

t c (t) = t + T f + T b + r(t + T f ) 

Future information 

t(t c ) = 9(t c -T b ) -T f 



(27) 



Past information 
RTT{<} = t c -t 

= tc-g(tc-T b )+T f 
= T b + T f +T(g(t c ~T b )) 

where Tf and T b are the forward and backward propagation 
delays corresponding to the constant delay operators Vo,i and 
Vi, 2 . 

Using the backward expression of RTT and the relation 

o £g Ci = S$q ° &Ci = Id, it easy to obtain the following 
result: 

Result 19. The flight size obeys 



^ (t) 

Fi{C u &cM = [<f>(uf,s)ds 
t 

t 

Fi{Ci,t) = J (f>(ut,s)ds. 
m Ci (t) 



(28) 



(29) 



Proof: This is an immediate consequence of the 
conservation-law (01 (through the ACK-clocking model) and 
the above RTT models. ■ 
Note that in ll22l . a model is obtained directly from the W- 
ACK-clocking applied directly to a selected topology. How- 
ever, the model is not modular itself since hand calculations are 
needed to make it implementable that results unfortunately in 
complexity very sensitive to the topology. Hence the objective 
of obtaining a metamodel is not attained. The reason is that the 
ACK-clocking model is used in implicit form while our pro- 
posed method incorporates solutions of it, yielding an explicit 
formulation achieving the characteristics of a metamodel, i.e. 
modularity and scalability. 



D. ACK-Clocking Dynamics and User Flow Computation 

Since the flight-size expression ( |29l is exactly of the form 
(fJJ, then Proposition [3] can be immediately applied to derive 
an explicit expression for the output flow of a given circuit, 
that is the flow of received acknowledgments: 

Result 20. Let us consider a circuit C{ — {uf , u~f ). Then, the 
ACK-flow that user m receives is given by 

4>(u-,t) = <% Ci (t)<j>i(ut, M Ci (t))- (30) 



Proof: The key idea is to remark that Fi(Ci,t) = 
N + (t, S§Ci{t))- Hence, using Proposition [3] and noting that 
the ACK-flow corresponds to the flow leaving the circuit 
4>(u~,t), the result is obtained. Differentiability of &d(t) 
is inferred from the differentiability of the backward delay 
operators. ■ 
Note that differentiation of ( 1291 yields 

<p i (ut,t) = F' i (C i ,t)+M Ci (t)(f> i (uf,g!c i (t)) (3D 

meaning that to maintain a constant flight size, i.e. F'i(Ci,t) = 
0, the user has to naturally send data at the same rate it receives 
ACK packets: this is exactly ACK-clocking but expressed at 
a flow level. By flow integration, we can easily recover the 
'packet-level ACK-clocking'. A similar expression is reported 
in IfTSl . but stated without proof and using the model in 
([Tol l to represent the buffer output flows. Once again the 
proposed methodology allows to provide strong mathematical 
foundations for some existing results. 

E. User Flow, Flight-Size and Congestion Window Size 

We need to clarify the relation between a user congestion 
window size Wi(t) and its sending rate <p(uf,t). First, recall 
that the congestion window size corresponds to the desired 
flight-size, while the flight-size is the current number of 
packets in transit. The window size is then a reference to track 
while the flight size is the controlled output. The control input 
is the user sending rate. 

When the window size increases, the user can immediately 
send a burst of packets to equalize the flight- and window- 
sizes. In such a case, we can ideally assimilate them to 
be equaj^l (and so are their derivatives). The small delay 
corresponding to the protocol reaction time can be easily 
incorporated in the constant part of the RTT. The problem is, 
however, slightly more difficult when the congestion window 
size becomes smaller than the flight-size. In such a case, we 
can not withdraw packets from the network and the only thing 
we can do is to wait for the packets in the network to be 
acknowledged until, at some point, the flight size becomes 
equal to the window size. This basically means that the 
rate of decrease of the flight-size is equal to the rate of 
received acknowledgments (the rate at which data leave the 
network). Therefore, while positive slope of the flight-size is 
ideally unconstrained from above, the negative slope is lower 
bounded. In fl9), a rate-limiter is a posteriori added to the 



7 This is the main assumption in 1221 justifying the use of the W-ACK-clock- 
ing model. 
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model in order to constrain the negative slope of the flight size. 
This solution is however difficult to implement in the context 
of |9] due to the time-varying nature of the slope lower bound 
and the lack of any ACK-fiow model. Note also that in most 
recent works lTT8l . ||221 . this problem is automatically excluded 
by considering that the flight-size is always smaller than the 
window size, and that the window size does not decrease 
'too-much'. To the authors' best knowledge no well-rounded 
solution has been provided yet for the problem of congestion 
window size decrease. We provide below an explicit and 
complete solution to this problem, regardless of the rate of the 
variation of congestion window size. This is achieved through 
an augmentation of the user model and the consideration of 
the flow of ACK packets. 



K*) - 

defined in OO 



User (u ,u + ) 



defined in d33b 



Fig. 6. User block 



According to the above discussion, the flight-size must obey 

«*(*) if Ti(t) 



F'iiGut) 



{u i ,t) otherwise 



where %{t) is a condition which is true when no lower limit 
on the rate of variation of the flight-size is imposed and false 
otherwise. 

Result 21. The flight- size Fi{Ci,t) satisfies ( 1321 ) if the user 
sending rate is defined as 



Wi{t) + 4>(uT ,t) if Ti(t) 







where %{t) = ([^(t) = 0] A [w t (t) 

Mt) = { °wi(t) + 4>{ur ,t) 
Moreover, this model is the simplest one. 



otherwise 



(33) 



(u i , t) > 0]) and 



if Ti(t) 
otherwise. 



(34) 



Proof: The ACK-buffer tt,, taking nonpositive values, 
measures the number of ACK packets to retain in order to 
balance the flight- and window-sizes. When the virtual buffer 
has negative state, i.e. 7r,(i) < 0, the arriving ACK-packets 
have to be retained until the state reaches 0. Once zero is 
reached, the user can start sending again until the window size 
decreases too fast, i.e. Wi(t) < ~<fi(u^ ,t). Substitution of the 
user sending rate defined by (1331 and (l34l in (t3Tb yields the 
flight-size behavior (f32t . To see that the model is minimal, it is 
enough to remark that both conditions in %{t) are necessary. 

■ 

In order to characterize the ACK-retaining mode, the ACK- 
buffer (l34l has to be adjoined to the protocol model (|2TT >. 
resulting in an augmentation of the state of the user model. 
The protocol behavior depends on the measurements pi(x t ) 
which are functions of the overall network state x t ; this state 
is discussed in more detail in Section I Villi 

We are now in a position to define user operators from (1331 . 



Definition 22. The user operatorUi(wi) : M+ — > M+ mapping 
the ACK-fiow 4>(u^ ,t) to the sending flow (j)(uf 7 t) is given 
by 

<t>(u + ,t)=U{w)<j)(u-,t) (35) 
where U(w) — diagj {Ui(wi)} and Ui is given in ( 1271 ). 

VIII. General Network model 

Modular and independent models for transmission channels, 
buffers and users have been developed in Sections [V] [VI] and 
IVHI respectively. In this section we summarize the obtained 
results in a compact form involving dynamical systems and 
operators, and properties of the model are discussed. Notably, 
correspondence of the proposed model with existing ones is 
emphasized/recalled. 

A. General model 

The general network model takes the form 



k(t) = Af(xt,<l>(u-,t),<Kb-,t)) 
Fi(C h t) = f So {t) <j)(uf , s)ds 



(36) 



with 



(32) 









U{w) 


0" 
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m + ■ 








Tlub 


























5(t) (37) 



where = col((/>(u + , t), </>(&+, t), 4>{u~ , t), <j>(b~ , t)) and 

x = co1(t, z, 7r) are the flows and the state of the network, 
respectively. The hybrid models for user and queue dynamics 
are described by the nonlinear (discontinuous) functional Af 
obtained from equations (O, (0, (fJTJ and d34"i >. The notation 
xt is here to emphasize that the evolution of the network 
state depends on past state values Il23l . Note that adjoining 
the flight-size expression is needed to obtain a finite number 
of equilibrium points. Indeed, since the user flow is computed 
from the derivative of the flight-size, the equilibrium infor- 
mation is lost and can only be recovered from the original 
expression of the flight-size. At equilibrium we indeed have 
F* = w* = RTT* 4>\ where RTT* and <p* are equilibrium 
values for RTT and the sending flow of user Ui, respectively. 

This model thus takes the form of a descriptor nonlinear 
hybrid positive time-delay system with state-dependent and 
constant delays about which many theoretical questions are 
open: well-posedness, existence of solutions, uniqueness of 
solutions, stability of solutions, etc. Note also that in this 
paper, we have not discussed about delay-derivative constraints 
whose violation may lead to severe well-posedness problems 
[24], such as nonuniqueness of solutions, existence of small- 
solutions, stopping solutions, etc. Some simple topologies 
have been considered in |[T7l where it is shown that delay- 
derivative may exceed one under certain conditions. For the 
moment, it is unclear whether for arbitrary topologies and 
under certain reasonable conditions, the delays perceived by 
the users always have derivatives smaller than one. This 
property is very suitable for analysis since many theoretical 
tools can only been applied when this condition is verified, 
e.g. Lyapunov-Krasovskii functionals |25| or certain integral 
quadratic constraints [26|. 
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B. Model Approximations 

The proposed framework includes explicit and seemingly 
exact expressions for every quantity of interest, but this was 
not historically the case. The sending rate model has always 
been a missing link in past formulations where ad-hoc models 
were considered. It has been recently shown in [9|, [22 1 that 
these flow models are actually approximation of the W-ACK- 
clocking model, which is turn an approximation of the FS- 
ACK clocking model considered in this paper. We summarize 
these remarks below for completeness. 

1) Ratio flow model: By making the approximation 
Fi{t) ~ Wi(t) in the FS-ACK-clocking model to get the W- 
ACK clocking model, we obtain 



(38) 



Wi(t) = / cf>(uj , s)ds 
Jm Ci {t) 



and using the right-square rectangle rule we get the following 
expression for the sending rates 

Wi(t) 



(39) 



T. + rig^t)) 

which is very similar to the ratio flow model for instance 
considered in HTT1, 171. QQ. 

2) Joint flow model: The joint flow model reuses the W- 
ACK-clocking model and by making a first order Taylor 
expansion on the implicit expression 

Wi(t)- / (p(uf,s)ds = (40) 



we obtain 

w t (t) + (T t + T(t + T/)){w t {t)-0(u+,t))^O (41) 



or equivalently 



Wi(t) 



T L + r(t + T. 



A 



i(t) 



(42) 



which is exactly the joint flow model considered in |9|, |19|, 
[21 1, [20|. Neglecting the derivative term yields the usual ratio 
model HI], 0, ED. 

3) Static model: The static link model assumes that sending 
rates are proportional to the derivative of congestion window 
sizes, making the relation between these sizes and queuing 
delays static. This model can be obtained by further approxi- 
mating the ratio model or using linearization and a (0,0) Pade 
approximation. In Section IIX-BI a more general proof for the 
static-link model is provided and it suggests that the static 
model has a much wider domain of validity, as experimentally 
emphasized in Ii271 . 

IX. The Single-Buffer/Multiple-User Topology 
with Delay-based Protocols 

The purpose of this section is two-fold: exemplify the 
modeling technique on a simple topology and prove that when 
some conditions on the topology are met, the proposed model 
reduces to a model involving a static-link model [27 1 . The 
proposed model hence allows to clarify the status of the static- 
link model ll27l by providing, for the first time, a mathematical 
proof for its domain of validity. 



To this aim, a single-buffer/multiple-user topology con- 
nected by lossless transmission channels is considered. The 
forward and backward propagation delays of user Ui are 
denoted by t( and T\ respectively. We propose to use the 
following generic model of any delay-based congestion control 
protocol as the user model 



where z t , w,(t), Ti = T/ + Tf and g l (t) = g(t - T*) are the 
state of the protocol, the congestion window size, the prop- 
agation delay and the backward delay operator respectively. 
The functions V% and Wi are defined as in (1211 1. 



(43) 



A. Multiple-User/Single-Buffer Topology Model 
The general model is given by (1361 with (l43l and 



f(t) 
a(t) 



c _ V(t) +5(t) - : 



Wi(t) + (j)(u~,t) 

Sl Ht) _ T f4>{ut,e)de 

/ m{t) + 4>{u~ ,t) 
\ ° 

[ c <k u +, 9 1 M-t/) 
Eti ^f, t-Tf) 



if C(t) 
otherwise 
if %{t) 
otherwise. 



if Ti(t) 
otherwise 



if C(t) 
otherwise 



(44) 

where 6(t) denotes the normalized cross-traffic 5(t) £ [0, 1). 
The topology of this model is completely described 
by (1371 with the operators B, U = diag^,{Wi}, 
T^ub = [V T 6 ... V T bJ , K bb = and K bu = 
V T / . . . V T / where is the constant delay operator 
with delay d > 0. 



B. Homogeneous Delays and No-Cross Traffic - The Static- 
Link Model 



As stated in Section IVIII-B31 the static-flow model can be 
obtained using various approximations, which suggest that the 
static-model is only valid locally. This however contradicts the 
results reported in [27] where it is emphasized that the static 
model may yield quite precise results over a wide domain. 
Note also that the static-link model has been invalidated in 
many scenarios, notably some involving very heterogeneous 
delays or cross-traffic, see e.g. [9 1 and examples of Section [X] 
In the following, we show that the static model can be exact 
when some conditions on the network topology are met. 

Assuming that the propagation delays are homogeneous, i.e. 
T? = T*, Tf = T\ i = 1, . . . , N, and that there is no cross- 
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traffic, i.e. 5 = 0. 

m 

*<(«) 
Fi{t) 

H u 7 > *) = 

a(t) 
9b(t) 



the model in (l44l) reduces to 

■-^(t) - 1 




<7(i-r b ). 



9i,(t)-T/) 

,t-Tf) ' 



if C(t) 
otherwise 
if 7J(i) 
otherwise. 

if 7I(t) 
otherwise 

if C(t) 

otherwise 



(45) 

Assuming further that the buffer is always congested (i.e. C(i) 
holds true for all t) and all users are active (i.e. the 7l(t)'s are 
all true) we obtain 



r(t) 



i{t~T f ) 



(46) 



after the substitution of sending flows in queue dynamics. 
Integrating the above equation from to t we obtain 



-(i) = c- 



i(t-T f ) -T 



(47) 



where we assumed r(0) = and Wj(0) = 0, i = 1, . . . , N. 
The additional constant term — T can be determined such 
that the above equation satisfies the equilibrium equation 
J2i w*(T + r*) = c. Equation (|47T i is exactly the static model, 
showing then its exactness for the single-buffer/multiple-user 
topology with homogeneous propagation delays and no cross- 
traffic. This result might be generalizable to the case of 
chained buffers and multiple users. The exactness of the static 
model over more complex topologies is an open question. The 
case of single-buffer with constant cross-traffic may also be 
analyzable. 

Thus, according to the model proposed in (l44l . the static- 
link model is exact in the single-link topology whenever 

• the buffers are permanently congested, i.e. C(t) holds true 
for all t > 0; 

• the propagation delays are homogeneous, i.e. t( = , 
T\ = T b , i = 1,...,N; 

> the cross-traffic is absent, i.e. 5 = 0; 

• the users are not in ACK-retaining mode, i.e. %(t) holds 
true for all i = 1, . . . , N. 

Compared to the justification of this model in 11221 . the 
proof developed above is much more insightful since no model 
approximation is made, only assumptions on the network 
topology. This shows that the static-link model has an ap- 
plication domain which is much wider than the ratio-link and 
the joint-link models, when the conditions on the topology 
are met. It is indeed valid in the nonlinear setting and it does 
not result from any approximation, just assumptions on the 
topology. 

This also shows that the proposed metamodel is able to 
provide theoretical justifications of a simpler model. This 
supplies a way for deriving proofs for validity domain of 
models. 



C. Homogeneous Delays and Cross-Traffic 

When cross-traffic is added to the problem, the overall pic- 
ture changes. The cross-traffic acts as a bandwidth limiter both 
in the networking and control terminology. Indeed, a nonzero 
5(t) reduces the maximal output capacity c, creating then sort 
of 'varying-output-capacity' c(l — S(t)) which reduces the 
bandwidth perceived by users. 

Result 23. In the congested mode, the queue model writes 

m 



cS (9bf(t)) 

c5{g bS {t)) + T^M l9b} {t) - Tf) 



(48) 



where gbf(t) = g(t — T b — T*) and cp(t) is the output flow 
at time t corresponding to the cross-traffic. The term <p(f) is 
responsible for the bandwidth reduction. 

Proof: Since the buffer is always congested and the users 
are not in ACK-retaining mode, we have 



and 



<pi(uf,t) = Wi(t) +4'i{u l ,t) 
' o5{9b{t))+Y, j <f>M, 9b {t)-Tfy 



(49) 



(50) 



Substituting the above expressions in the queue model and 
noting that J2j 0j( u 7'* ~ = c y^ e ^ s ^ resu lt- 



X. Model Validation 

We consider the scenarios given in J9), ll22l to validate 
the proposed model. The results obtained via NS-2 have 
been slightly shifted in time so that the congestion window 
variation times match. Unlike 11221 . the results are not shifted 
in amplitude and this causes small discrepancies. If, however, 
the NS-2 results were shifted vertically so that they match 
initial equilibrium values, then the curves would match almost 
perfectly. 

A. Single-Buffer/Multiple-Users 




Fig. 7. Topology of scenarios 1 & 2 

In this section, we consider the interconnection of two users 
through a single resource, as depicted in Fig.[7j The bottleneck 
has capacity c = lOOMb/s and the packet size including 
headers is 1590 bytes. The following scenarios from J9], 11221 
are considered: 

• Scenario 1: the congestion window sizes are initially 
Wi = 50 and u>2 = 550 packets, at 3s wi is increased to 
150 packets. The propagation delays are T\ = 3.2ms and 
T% = 117ms for users 1 and 2 respectively; see Fig. [8] 
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Scenario 2: the congestion window sizes are initially 
Wi = 210 and w® = 750 packets, at 5s W\ is increased 
to 300 packets. The propagation delays are T\ = 10ms 
and T2 = 90ms for users 1 and 2 respectively; see Fig. 

m 




Fig. 8. Scenario 1: Queue size Fig. 9. Scenario 2: Queue size 



We can see that our results fit well with the ones obtained 
by packet level (NS-2) simulations. Yet, this is not the case 
for the ratio-link, static-link, and the joint-link models when 
the transient-state is considered. The obtained results by the 
proposed model are identical to the results given in ||9l , 
E2l and obtained using the W-ACK-clocking model. This is 
expected since the W-ACK clocking is an approximation of 
the FS-ACK-clocking model defended in this paper. Note also 
that the approximation condition related to the window size 
increase is satisfied here. 



Fig. 11. Scenario 3: 
and queue 2 (bottom) 



queue 1 (top) Fig. 12. Scenario 4: queue 1 (top) 
and queue 2 (bottom) 



The proposed model is again able to capture the network 
behavior well and it retrieves the previous results reported in 
|9l , ll22l . This is again due to the equivalence between the 
ACK-clocking models in this case. 

We introduce now a constant cross-traffic x c i = c\/2 on the 
first link. Initiall>E we set tuj = 1200, w§ = 1600, tog = 5 
and we consider the following scenarios: 

• Scenario 5: The congestion window W\ is increased by 

200 packets at 10s; see Fig. [13] 
> Scenario 6: The congestion window W2 is increased by 

200 packets at 10s; see Fig. [TJ] 




B. Multiple-Buffers/Multiple-Users 



Fig. 13. Scenario 5: queue 1 (top) Fig. 14. Scenario 6: queue 1 (top) 
and queue 2 (bottom) and queue 2 (bottom) 




s+ 5- 

Fig. 10. Topology of scenarios 3 to 6 (<5 represents cross-traffic I/O nodes). 



Here, we consider the case of two buffers interconnected 
in series (see Fig. ITDb with capacities c\ = 72Mb/s and C2 = 
180Mb/s. The packet size including headers is 1448 bytes. 
The link propagation delays are 20ms for link 1 and 40ms 
for link 2. The total round-trip propagation delays are T\ = 
120ms, T2 = 80ms and T3 = 40ms for sources 1, 2 and 3 
respectively. Initially, the congestion window sizes are w\ = 
1600 packets, w® = 1200 packets and w® = 5 packets. The 
following scenarios from [9], [22] are considered: 

> Scenario 3: No cross-traffic and the congestion window 
wi is increased by 200 packets at 10s; see Fig. Q~T] 

> Scenario 4: No cross-traffic and the congestion window 
u>2 is increased by 200 packets at 10s; see Fig. [12] 



The obtained results are identical to the results obtained by 
the the NS-2 simulations. Notice the reaction time between the 
moment at which the congestion window size is increased and 
the moment at which the second queue sees the flow variation. 
This illustrates that the model captures well the communi- 
cation path and the order of elements (spatial and temporal 
topology). These characteristics are not directly visible in the 
simulation results given in 11221 since the curve steps seem to 
have been aligned on the same temporal cursor. 

Despite of results equivalence, the proposed metamodel 
enjoys interesting properties such as modularity and scala- 
bility, that the model reported in 11221 lacks. This is an 
important improvement over previous models that were not 
able to cumulate accuracy, scalability, modularity and other 
interesting properties. This will be discussed in more detail in 
the next Section. 



C. Decreasing the Congestion Window Size 

The models proposed in [18], [22] do not capture sudden 
decreases in the congestion window size that would cause 



8 This scenario is actually identical to the one in 1221 Section III.B.3], the 
initial values for congestion window sizes given in 1221 are incorrect. 
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the buffer to empty or become smaller than the actual flight- 
size, that is, smaller than the number of packets in flight. The 
proposed model does capture these phenomena since 1) the 
FS-ACK-clocking model derived from the conservation law 
involves the flight-size rather than the congestion window size, 
unlike in [22|; and 2) the user model implements an ACK- 
buffer to count the number of packets to remove from the 
network before starting to send again. Note that the derivation 
of the user model including the ACK-buffer has been made 
possible due to the availability of an explicit expression for 
the flow of acknowledgments (f30b . This makes it computable 
through an explicit solution for the queuing delay and the 
buffer output flows Q3, |28). In Q, the deer easing of 
congestion window size is handled by adding a rate limiter 
to constrain the (negative) slope of the queue size. This rate 
limiter is however rather difficult to characterize due to the 
time-varying nature of the lower-bound on the slope which 
depends on the received rate of acknowledgment and the 
network state, the former being unfortunately unavailable in 
the framework of the thesis J9]. 

Let us consider the single-user/single-buffer case where 
the total propagation delay is T = 150ms, the packet size 
including headers is 1040 bytes and the initial value of the 
congestion window size is w° = 500. A t = 5 seconds, the 
congestion window size is halved. We consider the following 
scenarios 

• Scenario 7: c = 12.5Mb/s and no cross-traffic; see Fig. 

El 

« Scenario 8: c = 25Mb/s and half capacity used by cross- 
traffic; see Fig. [161 




Fig. 15. Scenario 7: Queue size (top) Fig. 16. Scenario 8: Queue size (top) 
ACK buffer (bottom) ACK buffer (bottom) 

We can see that we obtain exactly the same results as NS- 
2 simulations (and the rate-limiter model reported in |l9"1). 
As desired, the ACK-buffer measures (counts) the number of 
packets to remove before starting to send again. 

XI. Related Work 

Despite of being unorthodox, it was to the authors' point 
of view to discuss the related works after presenting the 
proposed metamodel. Many congestion models have been 
proposed in the literature with different formalisms: time- 
domains, protocols, network topologies and other hypotheses. 
It is thus quite difficult to compare them directly and give a 
thorough discussion. Yet, we have identified some important 
characteristics which are summarized in Tables [i] |TT] and Hill 
Table U collects global characteristics of a network model 



such as the time-domain property to indicate whether the 
model is in discrete (DT) or continuous (CT) time, or the 
modularity property to indicate whether it is modular or not. 
With delay we indicate that the model considers it either as 
constant (Cst), time-varying (TV), state-dependent (SD), or 
asynchronous (Async) discrete-time system whose asynchrony 
(incorporating delays) depends on the state of the system. 
Spatial and temporal topologies denote whether the model 
share the same spatial and time topologies with the actual 
network, respectively. Generic indicates the ability of the 
model to describe the network for different protocols whereas 
rate model denotes from which formula the user sending rate 
is computed. The ratio flow model flj], El, 0, H3, ® 
is given by quotient of the congestion window size and the 
RTT The static flow model 11271 assumes that the flow is 
proportional to the derivative of the congestion window size 
and the joint flow model (gj, Q9), ED, |f20| assumes that 
the flow is given by the sum of the ratio-flow model and the 
congestion window size derivative. 

Tables [II] and [TTT1 examine the important subparts of the 
model, namely the user and buffer models, respectively. It is 
also important to mention that none of them has the structure 
of a metamodel having the precision of the considered one. 
In Table [ill homogeneous and heterogeneous delays denote 
how well a model captures the case of homogeneous or 
heterogeneous propagation delays, whereas cross and bursty 
traffic denote how well a model captures the presence of cross 
traffic and bursty traffic [29], respectively. Exact indicates that 
the transient-state behavior of a model matches well to the 
one observed in packet level simulations. Too fast and too 
slow mean that a model converges to steady-state faster or 
slower, respectively, compared to what is observed in packet 
level simulations. Faster states a relativity faster convergence, 
yet not as much as observed when it is too fast. In Table 
UTTl output flows indicates whether it is defined as aggregate, 
meaning all input flows are mixed, or split, where the buffer 
acts as a square MIMO system mapping a given input flow to a 
given output flow. FIFO characterization denotes whether the 
model captures the actual content of the queue and the FIFO 
behavior. With queuing delay model we indicate whether its 
effect on input flows are clearly defined. Frequency filtering 
effect denotes whether the model of buffer have a filtering 
effect on the input flows to the output flows. 

1 ) Tabled 1 In this model, while the delay is denoted as a function of 
time, it is implicitly defined from the state, without further investigation. An 
exact delay characterization as a function of the state has been obtained in [17]. 
2 The state-dependent delay is implicitly characterized by the ACK-clocking 
formula and the RTT expression. 3 The modularity of the model could have 
been developed but was not. 4 Modular but not able to represent any topology, 
essentially due to lack of any output flow model for buffers. 5 This modeling 
technique can be applied to any type of topology, however hand calculations, 
which result in poor scalability, are necessary. 6 The notion of flow is not 
defined in a discrete-time model. 7 The discrete-time model is derived from 
a continuous-time model using the (FS)ACK-clocking model. 8 The flow is 
implicitly defined as the solution of the (W)ACK-clocking model. 9 The flow 
is computed by solving the (FS)ACK-clocking model explicitly. 

2) Table [771 10 This model is the one introduced in (28) and also 
considered in the current paper. 11 In the single-buffer case the model is exact. 
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No 


No 


No 


No 


No 


No 


No 
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Rate model 


No 6 


Ratio 


Ratio 


Ratio 
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(FS)ACK 7 


loint 


Ratio 


No" 


(W)ACK 8 


(FS)ACK 9 



TABLE I 

Main characteristics of some existing models in the literature. 





User sending flow models 


Static 


Ratio 


Joint 


(W)ACK 


(FS)ACK+ 
ACK-Buffer 10 


Congestion window 


Characteristics 














Homogeneous delays 


Exact/Too fast 11 


Too slow 14 


Faster 12 


Exact 


Exact 


Increase 


Heterogeneous delays 
Cross -traffic 


Too fast 
Too fast 13 


Too slow 14 
Too slow 


Faster 12 
Faster 


Exact 
Exact 


Exact 
Exact 




Bursty-traffic 


Exact/Too fast 15 


Too slow 


Faster 16 


Exact 


Exact 




Homogeneous delays 


Too fast 17 


Too slow 14 


Too slow 14 


Too fast 17 


Exact 


Decrease 


Heterogeneous delays 
Cross -traffic 


Too fast 17 
Too fast 17 


Too slow 
Too slow 14 


Too slow 
Too slow 14 


Too fast 17 
Too fast 17 


Exact 
Exact 



TABLE II 

Comparison of different user flow models with respect to the accuracy of estimated queue sizes. 



Buffer descriptions 


Aggregate model {£} {9} 


Pseudo-queue model 1151. 1181 


Flow model 1 1 3 |, 16|, |17| 


Characteristics 








Output flows 


Aggregate 


Split 


Split 


FIFO characterization 


Lost 18 


No 19 


Yes 20 


Queuing delay model 


Lost 18 


No 21 


Yes 22 


Frequency filtering effect 


All-pass/undef. 23 


All-pass/Low-pass 24 


All-pass 25 



TABLE III 

Comparisons between different buffer descriptions. 



It is unclear whether this is also true for the multiple buffer case. In most cases, 
the model is too fast. 12 Accurate modeling of discontinuity/high slope of the 
queue when burst of data arrives at the queue input. 1,j The model is usually 
too fast. It might, however, be possible to refine it to account for constant 
cross-traffic 1 19|. 14 Unable to characterize the high slope of the queue due to 
the low-pass filtering effect of the resulting queue model, see Section IVLDl 
15 Exact in the case of homogeneous delays and single-buffer topology since 
the relation is static and it has an infinite bandwidth; see Section IdTBI Too 
fast otherwise. 16 Larger bandwidth than the ratio model but still limited since 
the bandwidth reduces as the queue size grows. 17 These models may withdraw 
packets from the network by injecting negative flows. 

3) Table \III\ 18 The information is lost since the output flows are 
aggregate. 19 The output flows are insensitive to content swapping. 20 The 
output flows are sensitive to any swap of information in the queue. 21 The 
exact queuing delay is actually difficult to assign to any output flow due to 
the structure of the model for the output flows involving queues. 22 Undefined 
due to the aggregate formalism. 23 All-pass when the sum of the input- flows 
does not exceed the maximal output capacity, otherwise the filtering effect is 
not really defined due to the aggregate formalism. 24 All-pass when the sum of 
the input-flows does not exceed the maximal output capacity, otherwise acts 
as a low-pass filter with bandwidth inversely proportional to the queue size. 
25 This is due to the direct-feedthrough structure of the model. The delay and 
nonlinear formula for the output flows only have a compressing/expanding 
effect on time and amplitude. 



XII. Conclusion 

This paper presents a modular fluid-flow network conges- 
tion control model to analyze communication networks with 
arbitrary topology. This modular metamodel is introduced by 
mathematically modeling network elements such as queues, 
users, and transmission channels, and network performance 
indicators such as sending/acknowledgement rates and de- 
lays. It is composed of building blocks that implement local 
mechanisms, some of which being ignored/unmodeled by 
some existing models in the literature. It is shown that the 
proposed model allows to recover existing models and brings 
a formal proof for their validity/invalidity as well as for their 
domain of validity. We present a novel classification of the 
previously proposed models in the literature and we show 
that the existing models are often not capable of capturing 
the transient behavior of the network precisely. Numerical 
results obtained from packet-level simulations demonstrate the 
accuracy of the proposed model. We plan to extend our work 
with modeling of data loss, such as packet drops, and the 
time-out mechanism at the user level. 
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Figure 5: Comparison with NS-2 simulations of the total number of packets 
counted at the inputs of the buffers described by the proposed flow-based model 
and the queue-based model. 
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Figure 8: Scenario 1: Queue size 
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Figure 9: Scenario 2: Queue size 
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Figure 11: Scenario 3: queue 1 (top) and queue 2 (bottom) 
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Figure 12: Scenario 4: queue 1 (top) and queue 2 (bottom) 




11 11.5 

Time [seel 



12 



12.5 



13 



0.058 

^0.056- 

o 

CD 

«L 0.054- 
" 0.052 

§ 0.05- 

^ 0.048 : 

0.046 _ 
9.5 



10 10.5 



11 11.5 

Time [sec] 



I I I I 


— Proposed Model 




NS-2 








1 


1 1 1 1 1 



12 



12.5 



13 



Figure 13: Scenario 5: queue 1 (top) and queue 2 (bottom) 
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Figure 14: Scenario 6: queue 1 (top) and queue 2 (bottom) 
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Figure 15: Scenario 7: Queue size (top) ACK buffer (bottom) 
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Figure 16: Scenario 8: Queue size (top) ACK buffer (bottom) 
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